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We theoretically study phonon-drag contribution to the thermoelectric power and hot-electron 
energy-loss rate in a Rashba spin-orbit coupled two-dimensional electron system (2DES) in the 
Bloch-Gruneisen (BG) regime. We assume that electrons interact with longitudinal acoustic phonons 
through deformation potential and with both longitudinal and transverse acoustic phonons through 
piezoelectric potential. Effect of the Rashba spin-orbit interaction on magnitude and temperature 
dependence of the phonon-drag thermoelectric power and hot-electron energy-loss rate are discussed. 
We numerically extract the exponent of temperature dependence of the phonon-drag thermopower 
and the energy-loss rate. We find the exponents are strongly suppressed due to the presence of the 
Rashba spin-orbit coupling. 

PACS numbers: 72.20.Pa,75.70.Tj, 73.21. Fg 
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I. INTRODUCTION 

There has been a rapidly growing interest on spin- 
orbit coupled low-dimensional electron systems like two- 
dimensional electron systems (2DES) formed at semicon- 
ductor heterostructure interface, quantum wires, quan- 
tum dotsi"— etc. The usefulness of the spin-orbit cou- 
pling in condensed matter system was realized after the 
proposal of spin field effect transistor by Datta and Das^ 
and thereafter, various interesting theoretical and exper- 
imental studies have been reported. One main goal is 
to control and manipulate the spin degree of freedom 
of charge carriers in nanostructures so that spin-based 
electronic devices^""— and quantum information process- 
ing technology can be developed. The coupling between 
electron's spin and orbital angular momenta naturally 
arises when one makes non-relativistic approximation to 
the relativistic Dirac equation. There are various types of 
SOI present in the semiconductor heterostructures. Most 
commonly used SOI is the Rashba SOI (RSOI)^ which 
is due to structural inversion asymmetry in semiconduc- 
tor heterostructures such as GaAs/AlGaAs. The RSOI 
is always present in the system intrinsically due to the 
internally generated crystal field. One can also increase 
strength of the spin-orbit coupling with the help of an 
external gate voltag o^'^" . The RSOI modifies various 
propertieaii"— of a 2DES including effective mass, uni- 
versality of spin polaron mass, polaron binding energy, 
mobility, etc. 

The electron-phonon interaction plays an important 
role in determining transport properties of a 2DES. There 
are mainly two kind of electron-phonon interaction in 
semiconductor heterostructures, namely deformation po- 
tential (DP) and piezoelectric (PE) scattering potential. 
The deformation potential is the change in potential en- 
ergy of an electron due to lattice deformation. In an 
inversion asymmetry host crystal, an electric polariza- 
tion is induced due to lattice vibration and the potential 
corresponds to electric polarization is known as piezo- 
electric potential. The electrons are scattered by the 



deformation and piezoelectric potentials and provides a 
non-zero contribution to momentum relaxation time, in 
addition to other contributions coming from disorders, 
impurities etc. Numerous studies have been devoted to 
probe electron-phonon interaction (EPI) by measuring 
low-temperature mobility^""— of a 2DES. 

Thermoelectric properties of various materials includ- 
ing 2DES have attracted much interest due to potential 
applications. With the application of an external tem- 
perature gradient VT across a sample, an electric field 
E (X VT is generated. The proportionality constant (5) 
between E and VT is known as thermoelectric power 
or the Seebeck coefficient. There are mainly two con- 
tributions to the thermoelectric power S: diffusion ther- 
mopower Sd and phonon-drag thermopower Sg. The ap- 
plied temperature gradient gives rise to flow of electrons 
and phonons from hotter region to cooler region. The dif- 
fusion thermopower is solely due to flow of electrons and 
sensitive to the energy dependence of various scattering 
mechanisms such as ionized impurity scattering, surface 
roughness scattering etc. On the other hand, flow of 
phonons will try to drag the electrons from hotter region 
to cooler region due to electron-phonon interaction and 
giving rise to phonon-drag thermopower. Extensive the- 
oretical and experimental investigations^"— of phonon- 
drag thermopower of a 2DES without RSOI have been 
performed. 

Another relevant mechanism for probing electron- 
phonon interaction is the energy-loss rate (P) of hot 
electrons. When an electron system is subjected to uni- 
form heating, electron temperature raises above that of 
the phonons. Hot electrons relaxes to lower temperature 
via acoustic phonon emission. There are numerous ex- 
perimental and theoretical studies'^-"^^ on hot-elecrons 
energy-loss rate of a 2DES without RSOI. 

Unlike diffusion thermopower, both phonon-drag ther- 
mopower and hot-electron energy-loss rate depend only 
the electron-phonon interaction. Therefore, one can de- 
termine the electron-phonon coupling constant reliably 
by measuring Sg and P. 
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Very recently, diffusive thermopower-^ and acoustic 
phonon-limited resistivity'*^ in a spin-orbit coupled 2DES 
have been studied. To the best of our knowledge, a de- 
tailed study of the effect of RSOI on phonon-drag ther- 
mopower and hot-electron energy-loss rate have not been 
studied yet. In this paper, we study phonon-drag contri- 
bution to the thermoelectric power and energy-loss rate 
of a Rashba spin-orbit coupled quasi-2DES in which two- 
dimensional electron wave vector k couples with three- 
dimensional phonon wave vector Q — (q, f^z). We con- 
sider both DP and PE scattering mechanisms responsible 
for the electron-phonon interaction. In the BG regime 
we find analytically that Sg is proportional to T'* and 
for DP and PE scattering, respectively. On the other 
hand, approximate analytical calculations show that the 
energy-loss rate is proportional to and for DP and 
PE scattering, respectively. However, our numerical re- 
sults reveal that the exponents are strongly dependent 
on the electron density and the Rashba spin-orbit cou- 
pling constant. 

This paper is organized as follows. In section II we 
present all the analytical results of phonon-drag ther- 
mopower and energy-loss rate. Numerical results and 
discussions have been reported in section III. We sum- 
marize our work in section IV. 



and 



V2 



(3) 



with A = ± indicates two spin-split energy branches and 
tan(/)/j = kx/ky. At a given Fermi energy ep, the Fermi 
wave vectors for the two energy branches can be writ- 
ten as kp = yj [k^pY ~ ~ with kp = %J2'nn^ and 
ka = m* ajt? . The velocity of an electron in a particular 
branch A is given by 



h dk 



hk 



(4) 



A. Phonon-drag thermopower 

We consider the interaction between electrons with 
two-dimensional wave vector k and acoustic phonon with 
three-dimensional wave vector Q. To calculate phonon- 
drag thermopower we follow the explicit formula given in 
Ref*^ for 2DES. With appropriate modifications the ex- 
pression for phonon-drag thermopower in a Rashba spin- 
orbit coupled 2DES is given by 



II. THEORETICAL DETAILS 

We consider a quasi-2DES formed at the interface 
of semiconductor heterostructures which has a finite 
thickness in the confining direction (say, z). Typi- 
cally, the confining potential in the z-direction is a 
triangular potential. We assume that only the low- 
est sub-band due to transverse confinement is occupied 
by the electrons. Therefore, electrons are restricted to 
move in the xy plane with wave vector k = {kx,ky). 
One can write the electron's wave function as V'(r) = 
'iJj{x,y)^Q{z). The Fang- Howard wave function^ in the 
z-direction is given by ^o(^) = \/b^ /2ze~''^/^ with b = 

(487rm*e^/eoK?i^)^'^(nd + llne/32j as the variational 

parameter. Here, m* is the effective mass of an electron, 
K is the dielectric constant. Also, £o is permittivity of 
free space, is the depletion charge density and rie is 
the density of electron. 

The single electron Hamiltonian is given by 
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-(To 



[C^xPy - CTyPx 



where p = ?ik is the momentum operator for the electron, 
ctq is the 2x2 identity matrix, a is the RSOI coupling 
constant and are the usual Pauli spin matrices. The 
energy eigenvalues and the normalized eigenstates cor- 
responding to the above Hamiltonian are, respectively, 
given by 



2j,2 



fi^k 
2m* 



Xak 



<^ 2aAkBT^ 



J2 E f^Q.n^t)h~fi4'') 



X' k,k',Q 



X I^^^'(k,k'){r(6^K-r(6^;R;}-v„ (5) 

where e is the electronic charge, Tp is the phonon mean 
free time, A is the area of the sample, a is the Drude 
conductivity, fcs is the Boltzmann constant, loq ~ VsQ, 
T(efe) is the energy-dependent momentum relaxation time 
of an electron, /(e) = [e^('^~t^'> -f 1]^* is the Fermi-Dirac 
distribution function with /3 — 1/(/cbT), is the ve- 
locity of an electron in a particular branch A, Vp is the 
phonon velocity defined as Vp = WsQ/Q and Wq^ (k, k') 
is the transition probability which is responsible for mak- 
ing transition of an electron from an initial state |k. A) to 
a final state |k'. A) with the absorption of a phonon. The 
explicit form of the transition probability is given by the 
Fermi's golden rule 



27r, 



W^' (k,k') = ^|C, 



iW' |2 



Nq5 



'{4 



'k',k+q, 



(6) 



(1) where {Cq" ^ is the matrix element responsible for the 
EPI and Nq = [exp{Phu!Q) — is the equilibrium 
phonon distribution function. 

With the help of Kronecker delta symbol (5k'.k+q in 
Eq. ([6]) the summation over k' in Eq. ([5]) can be easily 
evaluated by replacing k' with k -I- q. The argument of 
the delta function in Eq. ([6]) confirms the conservation 
of energy e^' = efe + hujQ ■ At this point we assume that 
T(efe) is approximately constant over an energy scale of 

(2) the order of hajq so that we can write T{£k+fiajQ) — T(efe). 
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The summation over k in Eq. ([5]) can be converted into 
an integral over by the following transformation 



E 



A 



(27r)2 h? 



Using the standard result / depe^NQ{NQ + 1) = 
(ri)!C(?T.)(fcBT)"+^ with C,{n) is the Riemann zeta func- 
tion, we obtain 



where 9 is the angle between k and q. At very low 
temperature we can make an additional approximation 
as /(efc){l - /(efe + ^wq)} ~ fuvqiNg + l)6{€k - ep). 
Now we convert the summation over Q into an integra- 
tion over q and Qz a.s J2q ~^ (I/^tt^) J qdqdqz. In BG 
regime, phonon energy is very small compared to the 
Fermi energy and consequently we can make a further 
approximation as q << 2kp. With all these assumptions 
described above taken into account one can obtain a fi- 
nal expression for the phonon-drag thermopower as (for 
intermediate steps see the Appendix A) 



where A — VsTp is the phonon mean free path. The total 
phonon-drag thermopower due to DP scattering is given 

by 



eTpVs'm*T(ep) 



gPP _ m*'AD^kB £^ 5!C(5)(fc,,r)4 ^^^^ 

The total phonon-drag thermopower Sg for DP scattering 
is proportional to T^. 

Calculations similar to DP scattering will yield the to- 
tal phonon-drag thermopower for longitudinal and trans- 
verse PE scatterings as 



dqdq^q' 



\Cc 



Q 



-ihlOQfNQ{NQ + l). (8) 



S. 



The matrix elements for DP and PE scatterings are 
respectively given by^ 



and 



PE _ 45m*2A(e/ii4)2fcB / e% 3!C(3)(fcBT)2 



(14) 



\C, 



_ D'^^Q l + AA'cos6i 



gPE _ 59v,tm*^A{ehufkB I e% i\a?,){kBTf 



Q \DP 



(15) 



and 



\C. 



X \' |2 (e/ii4)2?i 1 + AA'cosf 

X |/fe)pAi(t)(q,q^), 



3a, A' 



(10) 



where D is the DP coupling constant, hi4 is the rele- 
vant PE tensor component, pm is the mass density, Vsi(t) 
is the longitudinal (transverse) component of sound ve- 
locity, Aii(i,qz) = 9g2g4/[2(g2 + g2)3] and At(q,(j,) = 
(8(7fg2 ^ q<^)/[4{ql + q'^)^]. The Kronecker deha symbol 
(5a. A' in the matrix elements implies that the electron- 
phonon interaction is spin-independent. Finally the form 
factor 1/(52)12 -^^jjich is responsible for the finite thick- 
ness of the quasi-2DES and it has the form |/(<Z2)P = 
I / (iz^Q(z)e"?^^|2 — b^/iql + b'^Y foi" ^ triangular poten- 
tial. 

The phonon energy is given by = hujQ — 
hvg \Jq^ + q^, so we can write q — cos 4> / (hv s) and 
qz = ep sin (j) / {hvs) so that dqdqz — >■ Cpdcpdcj) / {hv s)"^ . We 
consider the quasi-2DES is very thin i.e. qz << b so 
|/((7z)|2 can be approximated as |/(gz)|2 ~ 1. With these 
substitutions Sg becomes 



Total phonon-drag thermopower due to PE scattering 
is given by S^^ = S^f + 2SPf and in this case S^^ - 
T2. 



B. Hot-electron Energy- loss rate 

The average energy-loss rate per electron via acoustic 
phonon emission is given by 



P 



dt) 



1 



A,Q 



ON, 



Q 



dt 



(16) 



where Ng is the total number of electrons. The rate of 
change of phonon occupation number for a given branch 
A is given by 



rdN, 



Q 



dt 



k 



W^\k, k + q){(7VQ + l)/(e^ + hwQ) 



S: 



eTpVsm*T{ep) 
^ ^ ~ Sn^h-iakBT^k^ 



1-A 



{hvs 



dtptlNQ^Nq + 1) [ d</>cos2 0. (11) 



X [l-/(e^)]-7VQ/(e^)[l-/(e^ + ^Q]}, 

(17) 

where the Fermi-Dirac distribution function /(e) is de- 
scribed by the temperature of hot electrons Tg and 
phonon distribution function is described by the lattice 
temperature T;. Obviously Tg is larger than T; so that 
electron can relax its energy via acoustic phonon emission 
and equilibrate to the lattice temperature. 
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Now we use the following identity 

1 - /(£k + hug) _ fjek + f^q) ^p^Hu^^^ 



1 - 



(18) 



with = l/iksTe). Using Eqs. ([TH) to ([18]) and taking 
all the assumptions made for calculating Sg into account 
we finally obtain the following expression for energy-loss 
rate as (detail calculations are given in Appendix A2) 



P = 



dqdq:^ 



with /3i — l/lkgTi). The integration over q and q^ in 
Eq. ([2]) can be easily evaluated by the same technique as 
described in the previous sub-section. The expressions 
for total energy relaxation rate for DP, longitudinal PE 
and transverse PE scattering are respectively given by 



put- _ 




{r: - Tf), (20) 



and 



PE _ 13m*2(e/ii4)2fc|2!C(3) 




-{T!-Tf) (21) 



(T^ - Tf ).(22) 



In the BG regime P^^ is proportional to for DP 
scattering and P^^ ~ for PE scattering. 



III. NUMERICAL RESULTS 

In the previous section we have presented approxi- 
mated analytical results of the phonon-drag thermopower 
and energy-loss rate. In this section we discuss numeri- 
cal results of phonon-drag thermopower and hot-electron 
energy-loss rate. To do this we solve Eqs. ([5]) and 
numerically for both DP and PE scattering mechanisms 
in the low temperature regime. For the numerical calcu- 
lations we consider material parameters of GaAs/AlGaAs 
heterostructures as m* = 0.067me with free electron 
5.31 X 10^ Kg m-3, vs = v,i = 5.12 x 10^ 



mass me, 

ms-\ Vst ^ 3.04 X 10^ ms"\ D = 12 eV, hu = 1-2 x 10^ 
,-1 _ioni ^._ini4_-2 = 10-11 eV m and 



Vm-\ K = 12.91, Ud = 10'^ r 
no — 10^^ m~^. 

We estimate the efi^ective exponent from the log-log 
plot of the phonon-drag thermopower versus temperature 
due to DP scattering in Fig. 1 for both a — and a = ao 
with different densities. At very low temperature (T ^ 1- 
3 K), we obtain v = 3.449,3.942 and 4.139 for Ue = 3no, 




T(K) 



T(K) 
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FIG. 1: (Color online) Plots of the phonon-drag thermopower 
due to DP scattering versus temperature for different values 
of the density. Here, solid, dotted and dashed lines represent 
Ue = 3no, Ue = 5no and rie = 7no, respectively. 




FIG. 2: (Color online) Plots of the phonon-drag thermopower 
due to PE scattering versus temperature for different values 
of the density. Here, solid, dotted and dashed lines represent 
Tie — 3no, rie — 5no and rie = 7no, respectively. 



5710 and 7no, respectively, when a — 0. On the other 
hand, we obtain ly = 2.876, 3.268 and 3.534 for Ug = Suq, 
5no and 7no, respectively, when a = ao. The value of v 
with a — ao gets lowered than that with a = 0. Figure 
1 also depicts that the magnitude of Sg with a = ao is 
less in comparison with a = 0. 

In Fig. 2, the log-log plot of the phonon-drag ther- 
mopower versus temperature due to PE scattering is pre- 
sented for different values of a and Ue- When a = 0, 
we obtain i' = 1.688, 1.961 and 2.096 for rig — 3no, 
5no and 7no, respectively. When a — ao, we obtain 
1/ = 1.290, 1.520 and 1.658 for rie = 3ri0: 5no and 7rio, 
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respectively. Comparing Fig. 1 and Fig. 2, one can con- 
clude that the magnitude of Sg due to DP scattering is 
larger than that due to PE scattering. 
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FIG. 5: (Color online) Plots of the energy-loss rate due to 
DP scattering as a function of electron temperature Te for 
a — and a = ao for different density. We set the lattice 
temperature Ti = 0. Here, solid, dotted and dashed lines 
represent rie = 3no , rie = 5no and rie = 7no , respectively. 



FIG. 3: (Color online) Plots of the phonon-drag thermopower 
due to DP and PE scattering versus temperature for a fixed 
density = 5no. Here, solid and dashed lines represent 
a = and a — ao, respectively. 



In Fig. 3 we plot Sg due to DP and PE scattering 
versus T for a = and a = ag at a fixed density rig = 
5no. Figure 3 clearly shows that slope of the line with 
a = ao is less than that with a = case. 




FIG. 4: (Color online) Plots of the phonon-drag thermopower 
due to DP and PE scattering versus a for different values of 
the density at fixed temperature T = 2K. Here, solid, dotted 
and dashed lines represent rie = 3no, rie = 5no and = 7no, 
respectively. 



Variation of Sg due to DP and PE scattering as a func- 
tion of the the Rashba coupling constant (a) is shown 
in Fig. 4. The phonon-drag thermopower Sg decreases 



rapidly with the increase of a in the case of DP scatter- 
ing. However, Sg due to PE scattering decreases with 
increase of a very slowly in comparison with the case of 
DP scattering. 

It is important to compare the orders of magnitude 
of diffusion and phonon-drag contributions to the total 
thermopower. The diffusion thermopower^^, is given by 

= (23) 
3|e|eF V cf^ 

where the parameter p depends on the scattering mech- 
anisms. We calculate Sd from Eq. (^5)) and Sg from 
numerical evaluation of Eq. ([5]). With a = ao and 

n = 5no at T = 2 K we obtain Sd 5.464 ^V/K, 

S^P 18.176 AiV/K and S^^ 7.135 /iV/K. It is 

clear that at T = 2 K the phonon-drag thermopower 
due to DP scattering dominates over both PE scattering 
induced phonon-drag thermopower and diffusion ther- 
mopower. 

Now we turn to present the numerical calculations for 
energy-loss rate of hot electrons in the BG regime. In 
Figs. 5 and 6 we have shown the dependence of P with 
electron temperature Te for DP and PE scattering, re- 
spectively. In general, energy-loss rate of hot-electrons is 
given by P = T{T^ - Tf) with F as the proportionality 
constant. By taking the lattice temperature T; = we 
determine the values of 6 for a = and a = ao from Figs. 
5 and 6. For DP scattering with a = the values of S 
are S = 4.707, 4.916 and 4.998 for rig = 3no, Buq and Tuq, 
respectively. When a = ao we obtain S — 4.268,4.511 
and 4.661 for Ue = 3no,5no and 7no, respectively. For 
PE scattering we find S = 2.945,3.078 and 3.139 for 
rie — 3no,5no and 7no, respectively at a = 0. With 
a = ao the values of 6 are obtained as S = 2.635, 2.790 
and 2.880 for = 3no, 5no and 7no, respectively. 

Our numerical calculations reveal that the effective ex- 
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FIG. 6: (Color online) Plots of the energy-loss rate due to 
PE scattering versus electron temperature Te for a = Q and 
a = ao for different density. Ifere, lattice temperature is fixed 
to Ti — 0. Solid, dotted and dashed lines represent = 3no, 
rie = 5no and rie = 7no, respectively. 




FIG. 8: (Color online) Energy-loss rate due to PE scattering 
is plotted as a function of Te — Ti with a = qq for two lattice 
temperatures Ti = 3.4 K and Ti = 5.5 K. Here, solid, dotted 
and dashed lines represent Ue = 3no , rie = 5no and Tie — 7no , 
respectively. 
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FIG. 7: (Color online) Energy-loss rate due to DP scattering 
is plotted as a function of — Ti with a — ao for two lattice 
temperatures Ti — 3.4 K and Ti = 5.5 K. Here, solid, dotted 
and dashed lines represent = 3no, = 5no and Tie = 7no, 
respectively. 
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FIG. 9: (Color online) Plots of the energy-loss rate due to DP 
and PE scattering versus a for different values of the lattice 
temperature Ti. We fix the electron temperature at Te = 6 
K. Here, solid, dotted and dashed lines represent Ti — 3.4 K, 
T; = 4.2 K and = 5.5 K, respectively. 



ponents of temperature dependence of phonon-drag ther- 
mopower and energy-loss rate strongly depend on the 
electron density and the Rashba spin-orbit coupling con- 
stant. 

In Figs. 7 and 8 we plot P due to DP and PE scat- 
tering as a function of Te — T; with a ~ ao for different 
density. Different values of T; have been considered here. 
Magnitude of P is higher at higher values of Ti and it de- 
creases with increase of the electron density. Comparing 
Figs. 7 and 8, we see that the energy-loss rate due to DP 
scattering is much higher than that due to PE case. 

In Fig. 9, we present energy-loss rate due to DP and 
PE scattering versus a for different lattice temperatures. 
It shows the energy-loss rates increase monotonically 
with a in both the cases. This behavior is quite differ- 
ent from the behavior of the phonon-drag thermopower 



versus a shown in Fig. 4. 



Density 


DP 


PE 


K) 


Q = 


a — ao 


Of = 


a = ao 


3no 


3.449 


2.876 


1.688 


1.290 


5no 


3.942 


3.268 


1.961 


1.520 


7no 


4.139 


3.534 


2.096 


1.658 



TABLE I: The effective exponent of the temperature depen- 
dence of Sg in the BG regime for various values of and 
a. 
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Density 


DP 


PE 


(ne) 


a = 


a — cto 


a = 


a — ao 


3no 


4.707 


4.268 


2.945 


2.635 


5no 


4.916 


4.511 


3.078 


2.790 


7no 


4.998 


4.661 


3.139 


2.880 



TABLE II: The effective exponent of the temperature depen- 
dence of P in the BG regime for various values of rig and 
a. 



IV. SUMMARY 

In this section we are summarizing the main results 
of the present work. In BG regime phonon-drag ther- 
mopower and hot-electron energy-loss rate have been 
calculated for quasi-2DES formed at the interface of 
GaAs/AlGaAs heterojunction. Both DP and PE scatter- 
ing mechanism have been taken into account separately. 
It is shown that the effective exponent of the temper- 
ature dependence of the phonon-drag thermopower and 
energy-loss rate strongly depend on the electron density 
and the Rashba spin-orbit coupling constant. For DP 
and PE scattering in the BG regime, the values of the 
effective exponents of T of thermoelectric power and hot- 
electron energy-loss rate for different values of Ue and a 
are summarized in Table I and Table II. 

It is shown that the orders of magnitude of phonon- 
drag thermopower due to DP and PE scattering are 
almost same and the phonon-drag thermopower domi- 
nates over the diffusion thermopower in the BG regime. 
The phonon-drag thermopower due to DP scattering de- 
creases fast with the increase of a. On the other hand, 
the phonon-drag thermopower due to PE scattering de- 
creases very very slowly with increase of a. 

The order of magnitude of the energy-loss rate due to 
DP scattering is much higher than that of the PE scat- 
tering. The energy-loss rate due to DP and PE scattering 
increases monotonically with a. 



Appendix A 

In this appendix, we briefly sketch the derivation of 
the phonon-drag thermopower. Using Eqs. ([5]) and (0 
we can write 



1 - A 



e/c + 



Q 

dOhwQW^^ (k, k + q)/(e^) { 1 - /(e^ -f /i^q) } 

{4){^t-^t+^)-^p- (Al) 



Using the approximation /(efc){l — /(efe -I- hug)} ~ 
hujQ{NQ + l)5(efc — ep), the integration over ek in Eq. 



(jAip can be easily done and we obtain 
_ eTpVsm*T{eF) 



1 - \. 

F ^ V 



X (iVQ + l)M/^\kp,k;^ + q)(v^^-vt+q) •Vp.(A2) 

Using Eq. dU and taking q « kp we can write the 
difference between two velocities as 



^ _ A- (— - ^^^^^ ' ^) 

m* hxkp kp 



.(A3) 



Similarly, the difference between two energies can be 
approximated as 



/fi kp , \ 

h Aa g cos a 

\ m* / 



m 
2m* 



where is the angle between k^^- and q. Converting the 
summation over Q into an integration over q, q^ and 
doing the 9 integration we can simplify Eq. (jA2l) as 



g\ ^ eTpV,m*^T{ep) / 
STT^h^akBT^k^^ ^ 

X (hUQfNQiNQ + l) 

a f m*Vs Q 
h\ hkp q 

with 

G(4,o) = {i- 



A 



dqdq^q^ 



QGik^p,Q) 



a 
+ T 
m" n 



(i-A,r^^)-4-v 

V \j£F + eaJ 2kV 



(A5) 



■m*Vs Q 



1-A 



— ) 



ep + e, 



With the assumptions q << 2kp and m*Vs << hkp we 
have G{kp,Q) ~ 1. Equation (fTT|) can be derived easily 
from Eq. (jASP with all the above mentioned approxima- 
tions taken into account. 



Appendix B 

In this appendix, we briefly we outline the derivation 
of the energy- loss rate of hot electrons. Using Eq. ^ 
and after doing the integration over 6, Eq. p7|) can be 
re-written as 



at 



nfr'q 



Ak 1- 



Ck + 



:NQiTi) 



k^G{k^,Q) 



(Bl) 
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At low temperature we use the same approximation 
made in Appendix A. After doing the integration over Ck 
in Eq. (jBip finally Eq. (fT6|) becomes 



1 - A 



\Cr 



Nq{Ti) 



f 



dqdqzhujQ 

(B2) 



By taking the approximation G{kp,Q) ~ 1 into ac- 
count at low temperature and after evaluating the inte- 
grations over q and qz in Eq. ()B2|) it is not difficult to 

get Eqs. (EnD-dlSD- 
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